Semi-classical Analysis of Spin Systems near Critical Energies 
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The spectral properties of su(2) Hamiltonians are studied for energies near the critical classical 
energy e c for which the corresponding classical dynamics presents hyperbolic points (HP). A general 
method leading to an algebraic relation for eigenvalues in the vicinity of e c is obtained in the 
thermodynamic limit, when the semi-classical parameter n _1 = (2s) -1 goes to zero (where s is the 
total spin of the system). Two applications of this method are given and compared with numerics. 
Matrix elements of observables, computed between states with energy near e c , are also computed 
and shown to be in agreement with the numerical results. 
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Introduction. Bohr- Sommerf eld (B-S) quantization 
formulae for nonregular values of the energy parameter 
have been set up in [1, 2] in the case of quantum Hamilto- 
nians in the Schrodinger representation. They differ from 
the regular case and show a logarithmic accumulation of 
the spectrum near energies corresponding to hyperbolic 
fixed points. In this letter we construct a general method 
to overcome the breakdown of standard B-S quantization 
near singular spectral points in the case of spins systems, 
a situation non covered by earlier results. 

su(2) Hamiltonians arise naturally in many areas 
of physics, in the study of mutually interacting spins 
1/2 systems or due to symmetries present in collective 
bosonic Hamiltonians. For such models the analog of the 
semi-classical limit is obtained when the number of in- 
teracting sub-systems increases (thermodynamic limit). 
The semi-classical parameter h is replaced by the the in- 
verse of number of interacting sub-systems n" 1 which is 
related to the size of the considered su(2) representa- 
tion, n — 2s. A typical example is given by the Lipkin- 
Mcshkov-Glick (LMG) model proposed in 1965 to de- 
scribe shape phase transition in nuclei [3]. This model 
is used to describe magnetic properties of molecules [4], 
interacting bosons in double-well structures [5] and to 
investigate the role of entanglement in quantum phase 
transitions (QPT) [6]. 

QPT [7] arising at zero temperature are related to 
non-analyticities of the ground state function of 

the Hamiltonian coupling constants. Since the non- 
analyticitics involved are generically algebraic, this kind 
of phase transitions are characterized by a set of critical 
exponents describing how physical quantities (density of 
states, excitation gap, observables) behave in the vicinity 
of such points. Recently, non-analyticities arising within 
the spectrum have received much interest [8, 9], they 
can be viewed as QPT arising for excited states [10]. In 
the semiclassical limit, this phenomena corresponds to a 
change in the topology of classical orbits and the appear- 
ance of hyperbolic points (HP). The non-analyticities in- 
volved are generically found to be logarithmic [8, 9, 10]. 



We derive and test numerically spectral analytical ex- 
pressions for two different situations arising as a LMG- 
like model (LMG plus a cubic term) where hyperbolic 
trajectories, homoclinic and heteroclinic, exist. Finally, 
we compute, both analytically and numerically, the ma- 
trix elements of observables between states near critical 
energies. We discuss their semi-classical behavior hoping 
to clarify critical phenomena arising at the so-called ex- 
cited states QPT. 

B-S Quantization and WKB. The non-normalized 
spin coherent states [11] for a su(2) representation of di- 



mension 2s + 1 are defined by \a) = e 



-s), a G C, 



with s, integer or half integer, being the total spin (in 
the following we set n = 2s). They form an over- 
complete basis with a resolution of the identity given 
by Jd/^gg = 1, where d„ = -Ms^M ^ and 
{a\a) = (l + aa) n . In the coherent states basis, the su(2) 
generators (S± = S x ± iS y ) act as differential operators 



na — a 2 d a 



d & ; S z = --+ad a , (1) 



on the space of polynomial functions iff (a.) = (a\iff) of 
degree n. To a generic operator 
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(2) 



where the PiS are polynomials in a, is associated a 
function (symbol) A(a, Q = 5Z,Pi(<5)C ■ I n the frame- 
work of the WKB approximation, eigenstates of an Her- 
mitian operator, iJ^(cv) = e^(a), are obtained set- 
ting ^(a) = e n ^ a o G ( Q ) dQ anc [ solving perturbatively the 
Riccati-like equation for G 

H [a,G(a) +n- 1 d & ] = e, (3) 

in powers of the semi-classical parameter n~ l , setting 



H 



i=0 



'Hi ; G = J2 n^Gi ; s = J2 n-%. (4) 
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The result is the WKB solution [12], given in terms of 
Hi [a,G (&)], 



*wkb(«) 
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x [1 + 0(0]. 
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Quantization of the energies is obtained by imposing that 
ty(a) is a univaluatcd function of a € C, implying that 
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(6) 



with fc e N, for all closed path 7. In the semi- 
classical limit the probability amplitude (a|a) _1 |\&(a)| 2 
of finding the system in the coherent state \a) is 
exponentially localized on the classical trajectory 

Co = {a : Ho (a, Yfkz) = £ o}, along which 
G 




FIG. 1: Phase space portrait of a classical trajectory Co (full 
lines) describing a critical orbit that passing through two HP 
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U. Moreover, if C contains no HP, the WKB Q * and "i- For °( n V2 ) < l Q 7 "*! < °( n °) both > the lin " 



solution (5) is an analytic function of a in its vicinity. In 
this case Eq. (6) can be explicitly computed by choosing 
7 = Co and using the semi-classical expansion of G. The 
result is the Bohr-Sommcrfeld quantization condition for 
a spin system [12, 13, 14]: To + n T\ + 0{nT 2 ) = 



where Tt 







-to L da-. 



= n k; 
L oj is the classical 
action obtained by integrated the symplectic 2-form 
u> = -U(l + aa) _2 d<5 A da over the interior of the classi- 

27T1 V ' 

e 1 -H 1 + ±d ( d 6l Ho 



cal trajectory E, and X\ = \ — 5^ 



27ri JC, 



da 



earized solutions around HP point and the WKB solutions, 
coexist (dark gray region), permitting to identify both asymp- 
totic behaviors. The "in" and "out" solutions are connected 
via the matrices. Branch cuts of the WKB solutions are 
displayed as broken lines. 



a second order differential equation, these four functions 
are obviously not independent. The explicit form of the 
Parabolic Cylindrical functions provides a "connection" 
between different asymptotic regions 



Quantization near HP. If an HP exits along 
classical trajectory, i.e. d(Ho = for some a, G Co, 
X\ diverges and the quantization condition has to be 
modified for energies of order n~ l around the critical 



energy eo — s c . Near such points, setting (3 
can be linearized and brought to the form 
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(8) 
(9) 
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with c = 



^2e_^ cosh^e-'M^+^r (ir? + 



by a simple transformation ^(a) = e np ^'^>(p), where 
p is a second order polynomial of j3. The constants Tfc 
depend on the parameter of the Hamiltonian around the 



HP. The solutions of 



*C9) 



are 



H(fi,nr 

given explicitly in the form of Parabolic Cylindrical func- 
tions [15]. Let us consider the following linear combina- 
tions of these two independent solutions, having a well 
defined behavior when \(3\ n 1 / 2 — > 00, for (3 in a vicinity 
of Co (see Fig. 1 for the directions along which each limit 
is taken), 



*out,_R(/3) 
*out,iC8) 

where p = 



/9" 



1 + O^l^n- 1 / 2 ) 

l + od^r^- 1 / 2 ) 



Til 

4- 2 



1/4 



and 77 = £ \ p -i T ™ ■ Being solutions of 



Constrains of the type (8) give a set of local rela- 
tions between the "in" and "out" basis. A set of non- 
local relations is obtained by identifying the asymptotics 
of WKB solutions (see Fig. 1), leading to "Fout.L = 
e 2mnS{an,a.j)^,.^ R ^ where S{a.i,a.j) is regularized action 
integral given in Table I, Vj = (±ir/j — h) depending on 
the side R/L and j indexing HP. ln(x) is defined as hav- 
ing a branch cut along the negative real axes, <7fc = 0, ±1: 
if the classical orbit does not cut the branch-cut of 
In (a — otk) and ±1 if it cuts it in the up-down or down-up 
directions respectively. Summarizing the local and non- 
local basis relations: \& ut = T\&i n and ^ ou t = rSPm; 
where V&out and 'i'm are column vectors collecting the 
"in" and "out" solutions for each HP (z), T and T are 
matrices, the first coupling states with the same (i) and 
the second coupling sates with (i) and (j) linked by the 
classical trajectory. Quantization is obtained by impos- 
ing the compatibility relation D = det(T — T) = 0. 



3 



2mS(a>i, ctj) = 27ri5 , 4j + ^(fj ln[(-l) CT J (a, - aj)] - v t ln[(-l) <Ts (a., — oii)] + a l mv l - Oj\m)A; 
o -a r&i a j - l r&i q r/- - - \/~* i M(— l) CT * («— M(— l)^ C^— "t)] j - 

27rig '.i = f Si T+^ da - - Jg, gg [(« - - q ^)Gi] - a< - fi ,. 1 dQ ; 

. 27ri5(Qi,ai) =27riSi + i^; 

2^ = /£ TTkd« - I ln[(-l) CT (a - [(a - S,)Gi] da; 

TABLE I: Regularized Action Integrals 



We now apply the general method presented above to a 
particular spin Hamiltonian 



H 



hS z 



7*^ 



(10) 



The Lipkin-Meshkov-Glick (LMG) model [3] is obtained 
from Eq. (10) setting = 0. The cubic term in Eq. (10) 
is added to provide asymmetric orbits in order to test the 
quantization relations in a case as generic as possible. For 
the LMG model a detailed analysis of the phase space and 
the characterization of the critical points can be found in 
[9, 16, 17]. For small values of \x the phase diagram pre- 
sented in [!)] is kept invariant. In particular the system 
conserves a homoclinic HP at a = for e c = — \h\ when 
7a; > \h\ < Ijyl and a heteroclinic caustic joining two HP 

for "f x > 7 y > \h\ corresponding to e c = 2 ^ y . For 

the homoclinic case, one obtains D = — cos ^ n ( SL + SR ^ _ 

sin { arg [r(l/2 - irfj\ + r) log ( 4p 2 ?i) + nn (S R - S L ) } , as 
in the case of the Schrodinger representation [2], where 
Sr/l are given by Si in Table I (directions of integration 
are given in Fig. 2). For the heteroclinic case the quanti- 
zation condition is rather lengthy and will be given else- 
where [ ]. The comparison of the semi-classical quan- 
tization conditions with numeric diagonalization of the 
Hamiltonian using a matricial representation of the spin 
operators is given in Fig. 2. In both cases the agree- 
ment between the numeric energies and the points where 
D = is remarkable, for the heteroclinic case one can see 
that the matching becomes worst as the modulus of the 
renormalized energy rj increases. 

Matrix Elements. In the semiclassical 

limit, the normalized matrix elements /^(e^ n ^) = 
(^m+fcm^m) Q £ an b ser vable A computed be- 

V / <* m +fcl* m + fc )<* m |* m ) 

tween eigenstates of an Hcrmitian operator H (with the 
energies £> m ) and e( m + fc ) ) ; are known to be simply given 
as the amplitude of the fc-th Fourier mode of the observ- 
able symbol A, evaluated along the classical orbit of en- 
ergyeM [20], f£(s^) = ± J% dt & k ^A[a(t), C(t)], 
where T is the period of the classical orbit and the flow 
equations are given by: d t a(t) = —id^H{a,Q; d t C,(t) = 
i <9a7i (a, Q) [Q is the variable conjugated to cv, for the spin 
case £ = ev(l + aa) -1 ]. This result stands for regular 
orbits and can be obtained considering the action-angle 
variables. Since / is the Fourier transform of a analytic 




FIG. 2: Homoclinic Case (Up): h — l;j x = 4; j y = 
1/4; fj, = 5. Heteroclinic Case (Down): h = l;7a, = 5; "/ y = 
2; fi — 6. Left: Comparison between the zeros of D (blue 
line) and the eigenvalues of (10) computed numerically (black 
dots) for n — 500. We define the renormalized energy 

h+e-, . „ _ 171+12 _ -(A+7 H ) v /7j7 



V = 



V = 



spectively for the homoclinic and heteroclinic cases. Middle: 
Stereographic projection of the critical classical orbit. Right: 
Critical orbit in Riemann Sphere, the zeroes of ^(a) (black 
dots) are plotted for n = 120, in the semiclassical limit they 
condense in branch cuts of Go [9, 19]. 



function the matrix elements vanish exponentially with 
increasing k. This is a generalization of the result early 
obtained by Heisenberg for the harmonic oscillator case. 

For singular orbits containing HP the period T di- 
verges, moreover no action-angle variables can be defined. 
Nevertheless it is still possible to estimate such matrix 
elements by analyzing local and global properties of the 
critical eigenstates [21]. Let us use the resolution of the 
identity in order to write matrix elements as integrals 
over Sj, a domain of size 0(n _1 Inn) around the HP <5;, 
and Sj j , a domain of order n" 1 around Cq. Within these 
two sets of domains the eigenstates are given, respec- 
tively, by special functions and WKB approximation, 



(^ m +k\A\^ m ) 



E / +E, 

(i)->Ci) (0 



r. 



The last equality follows from considering the symbol 
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A constant on the domain E,-, by orthogonality of the 
cigenstatcs this term is nonzero only for k = where 
it gives the norm of the eigenstate inside the domain, 
fi(i). The regular functions g^j(u)) = J^dt A{t)e ituJ 
are computed using the flow equations on the branch 
i — ► j. Since fi{i) oc Inn, we obtain at leading order, 



ft (e (m) ) 



E(i)"4(a»>Ci)M(») 



E(^ (j )g^[^(e (m+fc) -g (m) )] 
E(»)M(«) 



(ii) 

.(12) 



Diagonal matrix elements (mean values of observables) 
are thus given as a sum of ponderate weights of the dif- 
ferent HP and depend on local properties of eigcnstates 
near this points. On the contrary, non-diagonal elements 
are given by the global properties of the classical orbit. 
Since g A is analytic, the matrix elements will decay 
exponentially as the energy difference increases, however, 
near the critical energy the mean energy spacing is of 
order n(e^ m+k ^ — oa fc In - n, meaning that the 

exponential decay in k becomes slower with increasing 
n (see Fig. 3). For an observable with A vanishing at 
the HP the amplitude of all matrix elements vanishes as 
O(hrV), for fixed k (Fig. 3) [18]. This has a simple 
semi-classical explanation. In the critical case the vol- 
ume of the phase-space corresponding to an energy band 
of order n _1 around e c is <3(n _1 ) for regions of type Ej ; j 
and 0(n _1 Inn) for Ej. However, for A vanishing at the 
HP, the relevant regions to compute the matrix elements 
are Ej j which, by Heiscnbcrg inequalities, can carry 
only a finite number of states O(n ) and not the total 
O(lnn) eigcnstates. The only way of conciliating these 
two facts is to take a quantized observable described by 
an O(lnn) x O(lnn) matrix whose elements vanish in 
the classical limit. 

Conclusion. We have presented a method for comput- 
ing semi-classical spectra associated to any number of 
heteroclinic junctions. Not only the expected average 
spacing ~ In -1 n is observed, but an algebraic relation 
is derived for eigenvalues near the critical energy. The 
method is fully general and applies to any su(2) Hamil- 
tonian. In order to test it in full generality we have 
added a cubic term to the standard LMG model breaking 
the quadratic underlying symmetry. The agreement 
with numerics is remarkable, especially considering the 
fact that the formulas are algebraically quite heavy 
in the case of two hyperbolic fixed points linked by 
heteroclininc junctions. We have also computed the 
matrix elements of observables, and show that their 
semi-classical behavior is universal, and different from 
the one in the regular situation. Moreover we have given 
a physical argument for the logarithmic vanishing of 
these matrix elements in the classical limit. 
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FIG. 3: Matrix elements of the operator s z = — between two 
states near the critical energy. e^ c ' is chosen to be the energy 
closest to the critical classical energy e c . The agreement of 
the numerical data (dots) with the predictions of Eq. (12) 
(circles) gets better for n big. The logarithmic downward 
shift as n increases is due to the fact that ~ In -1 n. 



We are grateful to R. Mosseri, J. Vidal for fruit- 
ful and stimulating discussions and to P. Vicira 
for a careful reading of the manuscript. PR 
was partially supported by FCT and EU FEDER 
through POCTI and POCI, namely via Quant- 
Log POCI/MAT/55796/2004 Project of CLC-DM-IST, 
SQIG-IT and grant SFRH/BD/16182/2004/2ZB5. 



[1 

[2 

[3 
[4 
[5 
[6 
[7] 
[8] 

[9 
[10 

[11 

[12 

[13 
[14 

[15] 



J. N. L. Connor, Chem. Phys. Lett. 4, 419 (1969). 

Y. Colin de Verdiere and B. Parisse, Commun. Math. 

Phys. 205, 459 (1999). 

H. J. Lipkin, N. Meshkov, and A. J. Click, Nucl. Phys. 
62, 188 (1965). 

D. A. Garanin, X., M. Hidalgo, and E. M. Chudnovsky, 
Phys. Rev. B 57, 13639 (1998). 

V. V. Ulyanov and O. B. Zaslavskii, Phys. Rep. 216, 179 
(1992). 

J. Vidal, G. Palacios, and R. Mosseri, Phys. Rev. A 69, 
022107 (2004). 

S. Sachdev, Quantum Phase Transitions (Cambridge 
University Press, 1999). 

W. D. Heiss, F. G. Scholtz, and H. B. Geyer, J. Phys. A 
38, 1843 (2005). 

P. Ribeiro, J. Vidal, and R. Mosseri, Phys. Rev. Lett. 
99, 050402 (2007). 

M. Caprio, P. Cejnar, and F. Iachello, Ann. Phys. 323, 
1106 (2008). 

A. Perelomov, Generalized Coherent States and Their 
Applications (Springer, 1986). 

J. Kurchan, P. Leboeuf, and M. Saraceno, Phys. Rev. A 
40, 6800 (1989). 

R. Shankar, Phys. Rev. Lett. 45, 1088 (1980). 

A. Garg and M. Stone, Phys. Rev. Lett. 92, 010401 

(2004). 

I. Gradshteyn and I. Ryzhik, Table of Integrals, Series 
and Product (Academic Press, New- York, 1980). 



5 



[16] O. Castanos, R. Lopez-Pena, J. G. Hirsch, and E. Lopez- 
Moreno, Phys. Rev. B 74, 104118 (2006). 
[17] P. Ribeiro, J. Vidal, and R. Mosseri, arXiv:0805.4078. 
[18] T. Paul and P. Ribeiro (2008), in preparation. 
[19] S. Nonnenmacher and A. Voros, J. Phys. A 30, 295 



(1997). 

[20] T. Paul and A. Uribe, Ann. Inst. Henri Poincare (A) 59, 
357 (1993). 

[21] J. R. Cary and P. Rusu, Phys. Rev. A 47, 2496 (1993). 



